The energy-level statistics in the core of a vortex in a p-wave superconductor. 
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In the presence of strong disorder, the statistics of quasiparticle levels in the core of a vortex 
in a two-dimensional p-wave superconductor belongs to the universality class B corresponding to 
the ensemble of orthogonal matrices in odd dimensions. This novel universality class appears as a 
consequence of the 0(2) spin symmetry of p-wave pairing. It is preserved in the presence of random 
disorder, of electromagnetic vector potential, and of an admixture of the pairing of opposite chirality 
in the vortex core, but may be destroyed by spin-orbit coupling and by Zeeman splitting. 

The indications of p-wave superconductivity in Sr 2 Ru04 jl[ stimulated the study of exotic properties of p-wave 
^\ , superconductors. The order parameter in this compound is expected to be the same as in the A phase of 3 He, 
<i±(k) oc z{k x ± iky). The direction of the vector z of the triplet orientation is fixed by the anisotropy of Sr 2 Ru04 to 
. be perpendicular to the Ru-0 planes. Because of the strong anisotropy, one may consider a two-dimensional model as 
a starting approximation. In two dimensions, the p-wave superconducting gap does not have nodes, which resembles 
conventional superconductors. However, many differences appear in inhomogeneous setups including boundaries, 
vortices, and impurities. In particular, impurities generate bound states with circular currents |^|; similar subgap 
states appear at the boundary and at domain walls a single-quantum vortex posesses a zero-energy state of 
topological origin . 

It is this last property that motivates the present work. The zero-energy state will be shown to survive a certain 
class of disordered perturbations, such as a random potential (modeling the impurities) or a random electromagnetic 
vector potential. If such a disorder is strong, the quasiparticle levels in the vortex core mix and, at the energy scale 
much smaller than the superconducting gap, they may be described by a random-matrix ensemble. In the present 
y—i , paper we identify the corresponding ensemble with that of orthogonal matrices in odd number of dimensions (class B 
t— I ' in Cartan's classification) Q . Thus this example completes the list of universality classes corresponding to the eleven 
families of symmetric spaces (see Table 1). The three Wigner- Dyson universality classes correspond to A, Al, and 
All series (unitary, orthogonal, and symplectic, respectively) 0. Three more classes (AIII, BDI, and CII) appear 
in systems with massless Dirac fermions, as a consequence of the chiral symmetry ||. Finally, four more classes 
(C, D, CI, and Z?III) were shown to describe mesoscopic superconducting systems, depending on the presence of the 
spin-rotation and time-reversal symmetries Q . In the present paper we demonstrate that the remaining eleventh class 
B appears in p-wave superconductors under topologically-nontrivial (vortex-type) boundary conditions responsible 
Ch , for the zero-energy level. 

Let us start with briefly describing the properties of the ensemble of orthogonal matrices in odd dimensions (class 
B). The Lie algebra so(2N + 1) consists of real antisymmetric (2N + 1) x (2 AT + 1) matrices. Its dimension is 
dim[so(27V + 1)] = N(2N + 1) and its rank rk[so(2N + 1)] = N. If a matrix A belongs to so(2N + 1), the Hermitian 
matrix H — iA has one zero eigenvalue, and the remaining eigenvalues form pairs (o>i, — w.;), i — 1,...,N. At 
low energies we may, without loss of generality, assume the Gaussian probability distribution for the Hamiltonian 
dP(H) = cxp(—TrH^H/2v 2 )Y[dHij, where v is a large cut-off energy (in our problem, v is of the order of the 
superconducting gap). Then the joint probability distribution for the (positive) eigenvalues uJi is of the conventional 
form §: 

N 

dP{cJi} = | J{uJi}\ He-^^du;,, (1) 
where J{uJi\ is the Jacobian of the diagonalization of the Hamiltonian, 

N 

iJH}i=nk 2 -^nN Q - ( 2 ) 

i<j i—1 

At energies much less than the cut-off energy v, the correlations of the quasiparticle levels u>i are determined solely by 
the Jacobian J{uji} [the expression (Q) follows from the explicit form of the roots £( fe ) of the Lie algebra so(2N + 1) 
and from |J{tt>i}| = Y[k I Si^fc^h the values of a and (3 may also be easily found from dimension counting]. For 
so(2N + 1) the parameters of the level statistics are (3 = 2, a = 2. The values of a and f3 for the universality class 
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B coincide with those for class C, and only the presence or absence of the zero-energy level distinguishes the level 
distributions in the two classes. 

Next, we shall prove that a single-quantum vortex in a two-dimensional p-wave superconductor obeys the statistics 
of the so(2N + 1) ensemble, provided the symmetry of the Hamiltonian preserves the zero-energy level. Consider the 
Bogoliubov-de-Gennes Hamiltonian 
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where ^ a are the electron operators (a is the spin index), V^r) is the external potential of impurities, A(r) is the 
electromagnetic vector potential, A^ (r) and A y (r) are the coordinate-dependent components of the superconducting 
gap. [In the bulk, the preferred superconducting order is one of the two chiral components r/± — A x ± Ay, but in 
inhomogcncuos systems, such as a vortex core, an admixture of the opposite component is self-consistently generated 
fiof . We account for this effect by allowing the two independent order parameters A x and Aj,.] Star (*) denotes 
the symmetrized ordering of the gradients p M and the order parameters A^ [definition: A* B = (AB + BA)/2]. At 
infinity, the order parameters impose the vortex boundary conditions: 
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A y (r — > oo , <f>) = iA^e 
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where r and <fi are polar coordinates. Plus or minus signs in the exponent correspond to a positive or a negative single- 
quantum vortex. For an axially-symmetric vortex with the chirality of the order parameter non-self-consistently fixed 
(A y = iA x ), without the vector-potential A(r) and without disorder V(r), the low-lying eigenstates of the Hamiltonian 
(||) has been found by Kopnin and Salomaa [Q. The spectrum is 
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with loq ~ A 2 /ep. This result should be compared with the spectrum of the vortex core in a s-wave superconductor 
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The common feature of the spectra in the s-wave and p-wave cases is the symmetry about zero energy. If we 
interpret holes in the negative-energy levels as excitations with positive energies (and with the opposite spin), then 
this symmetry implies that the excitations are doubly degenerate in spin: to each spin-up excitation there corresponds 
a spin-down excitation at the same energy. For a s-wave vortex, this degeneracy is due to the full spin-rotation SU (2) 
symmetry. The p-wave Hamiltonian (0) has a reduced spin symmetry. Namely, it has the symmetry group 0(2) 
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This non-abelian group causes the two-fold degeneracy of all levels (except for the zero-energy level(s) where the 
symmetry 0(2) may mix the creation and annihilation operators for the same state). This symmetry is crucial for 
our discussion. Note that we have not included in the Hamiltonian neither the spin-orbit term (U50 • [tr X p]), nor 
the Zeeman splitting H(r) • <x. Either of these terms would break the spin symmetry 0(2), which would eventually 
result in a different universality class of the disordered system (type D with non-degenerate levels), if these terms are 
sufficiently strong. 

The difference between the s- and p-wave vortices is the zero-energy level in the p-wave case. It has been shown by 
Volovik that this level has a topological nature || . Indeed, suppose we gradually increase disorder in the Hamiltonian 
(H). The levels shift and mix, but the degeneracy of the levels remains the same as long as the symmetry 0(2) 
is preserved. The total number of levels remains odd, and therefore the zero-energy level cannot shift if the final 
Hamiltonian is a continuous deformation of the original one (without disorder), i.e. if the topological class of the 
boundary conditions (|J) remains the same. 

Now we proceed along the usual lines of the random-matrix approach. Let us take the most random distribution of 
Hamiltonians within a given symmetry class. The only symmetry of the Hamiltonian (^) is the spin symmetry 0(2). 
The time-reversal symmetry is already broken by the vortex and by the pairing, and therefore neither the vector 
potential A(x) nor local deformations of A^ can reduce the symmetry of the Hamiltonian. When projected onto 
spin-up excitations 7^ = /[u(r)^|(r) + u(r)\& |(r)]d 2 r, the Hamiltonian for the two-component vector (u, v) takes the 
form: 
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In an arbitrary orthonormal basis of electronic states, this Hamiltonian may be written as a matrix 



From the hermiticity of the Hamiltonian, it follows that h) = h. From the explicit form of the p-wave pairing, 
A = — A T (it is here that the p-wave structure of the pairing is important; for s-wave pairing we would have A = A T 
instead). These are the only restrictions on the Hamiltonian (||). If we define 

the restrictions on the Hamiltonian (|J) are equivalent to the condition that the rotated matrix HJqHUq 1 is real 
antisymmetric, i.e. it belongs to the Lie algebra so(M), where M is the dimension of the Hilbert space (the same 
rotation of the Hamiltonian was used in Ref. || to identify the D universality class). 

The last step in our argument is to note that, under the vortex boundary conditions, the dimension of the Hamil- 
tonian (^) is odd, not even (this may be difficult to visualize from the particle-hole representation (||), but easier from 
the rotated Hamiltonian UqHU^ 1 ). Thus for a single-quantum vortex, we identify the space of the Hamiltonians with 
so(2N + l). 

A simple consequence of this result is the level distribution (0) with a = (3 = 2. Besides the zero-energy level, 



this distribution is identical to that of class C realized in s-wave vortices [12 LJ. This allows us to use the trick of 
mapping onto free fermions to compute any correlation function of the density of states (DOS) ||. In particular, the 
average DOS is 

/ / sv 1 sin(27ru>/u>o) r/ s 

(p(w)) = v +8(w), (10) 

loq 2ttu> 

where uq is the average inter-level distance. 

One more lesson from our analysis is the difference in the universality classes of p-wave mesoscopic systems from 
their s-wave analogues. In particular, a p-wave system without a topological zero-energy level would belong to the 
universality class D (f3 = 2, a = 0), in contrast to the class C for its s-wave counterpart. A detectable physical 
consequence of such a difference is an increase of the average DOS near the Fermi energy (in class D) as opposed to 
a suppression of DOS near Fermi level in class C [in class B, such a suppression is compensated by a (5-function at 
zero] . 

In the present work we do not discuss the microscopic derivation of the level mixing. We assume "strong level 
mixing" which allows us to use the random-matrix approach. On the other hand, the symmetry of the p-wave pairing 
is known to provide certain "spectrum rigidity" suppressing shift and mixing of the low- lying levels by impurities [pj . 
Thus to drive the system into the regime of "strong level mixing" may require a stronger disorder than in a s-wave 
vortex. A microsopic picture of the crossover to the disordered regime in a s-wave case was developed in ^ ^] , and 
its extension to the p-wave case will be a subject of future studies. Furthermore, a disorder is known to suppress the 
p-wave superconductivity, and the possibility to reach the required level mixing before destroying superconductivity 
is not obvious. A coexistence of sufficiently strong disorder and p-wave superconductivity may possibly be achieved in 
alternative setups such as disordered normal-superconducting (p-wave) hybrid devices. Besides strong level mixing, 
the only requirements on the system to exhibit type B level statistics are the p-wave symmetry of the pairing and the 
topological zero-energy level. 

One more approximation made in our model is neglecting the Zeeman splitting. This is a good approximation 
for strong type II superconductors (with n = A/£ 3> 1). In Sr2Ru04 the experiments indicate k ~ 2.6 ]l5| ], which 
implies that the Zeeman level shift is of the order of the inter-level spacing u)q. However, in a clean vortex, this shift 
is approximately constant for the low-lying levels, due to the large coherence length (kp£ 1) and to a smooth 
magnetic-field profile H(r). It is likely that this property will hold even in the limit of strong disorder, then the 
overall level distribution will be simply shifted by a constant energy. 

An important observation related to the Zeeman splitting of the energy levels is the fractional spin 1/4 of the vortex 
in the ground state. Indeed, in the absence of the Zeeman field, the multi-particle ground state of the vortex is doubly 
degenerate, with the z-component of the spin in the two degenerate states differing by 1/2. The Zeeman field splits 



the ground state as if the vortex were a particle with the z-component of the spin equal half the spin of the electron. 
Thus we conclude that the vortex has the z-component of the spin S z = ±1/4. It would be interesting to understand 
possible physcial implications of this effect. 

Finally, for the convenince of the reader, we have gathered the information about all eleven symmetry classes in 
Table 1. This table is compiled from Refs. |^-^| and contains the dimensions and ranks of the symmetric spaces as 
well as the parameters a and j3 of the joint probability distributions of the energy levels (they may be computed 
solely from the ranks and dimensions by a simple power counting). The present work provides the example of the 
physical system belonging to the class B (the last line of the table). Besides, one more symmetry subclass has not 
been studied so far in the context of mesoscopics: the odd-iV subclass of DHL In the work of Altland and Zirnbauer 
, the even-A^ subclass of Dili is represented by a s-wave mesoscopic system with time-reversal symmetry and with 
broken spin symmetry. The novel odd-iV subclass of Dili may occur as a topological modification of that construction 
or in a p-wave superconducting system without time-reversal-symmetry breaking, provided a zero mode is required by 
topology. Finding a physically plausible mesoscopic realization of the odd-iV Dili subclass is an interesting problem 
in the framework of the symmetry classification of mesoscopic superconducting systems. 
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Wigner-Dyson 


AI [GOE] 


SU(N)/SO(N) 


(JV-l)(JV + 2)/2 


N — 1 
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All [GSE] 


SU{2N)/S P (N) 


(N - l)(2N + 1) 


N — 1 


4 




AIII [chGUE] 


SU(p + q)/S(U(p) X U(q)) 


2pq 


P 


2 


l + 2(q-p) 


Chiral ensembles [||, p < q 
(l ~ P) — number of 
zero modes 


BDI [chGOE] 


SO{p + q)/SO(p) X SO{q) 


pq 


P 


1 


q-p 


CII [chGSE] 


Sp(p + q)/Sp(p)xSp(q) 


■ipq 


P 


4 


3 + 4(g - p) 


C 


Sp(N) 


N(2N + 1) 


N 


2 
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Altland-Zirnbauer |j| 


D 


SO{2N) 


N(2N — 1) 


N 


2 





CI 


Sp(N)/U(N) 


N(N+l) 


N 


1 


1 


Dili 


SO{2N)/U(N) 


N(N — 1) 


m 


4 


1 (even N) 


5 (odd N) 


one zero mode 
(example unkonwn) 


B 


SO{2N+l) 


N(2N + 1) 


N 


2 


2 


one zero mode 
(present work 
andRcl @j 



TABLE 1. Symmetric spaces and universality classes of random-matrix ensembles. 



P.S. At the final stage of the preparation of this manuscript, the author has learned about the recent work of 
Bocquet, Serban, and Zirnbauer p6fl , where the type-B level statistics in p-wave vortices has been pointed out. 

The author thanks M. V. Feigel'man for suggesting this problem, for fruitful discussions and for helpful comments 
on the manuscript. Useful discussions with G. Blatter, V. Geshkenbein, D. Gorokhov, R. Heeb, and M. Zhitomirsky 
are greatfully acknowledged. The author thanks Swiss National Foundation for financial support. 
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